Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 1 February 2008 



Time-series Spectroscopy of Pulsating sdB Stars II: 
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ABSTRACT 

We present the analysis of time-resolved spectroscopy of the pulsating sdB star 
PG 1605+072. From our main observing run of 16 nights we have detected velocity 
variations at 5 frequencies that correspond to those found in photometry. Based on 
these data, there appears to be change in amplitude of the dominant modes over about 
a year. However, when we include extra observations to improve the frequency reso- 
lution, we find that some of the frequencies are split into two or three. Simulations 
suggest that the apparent amplitude variation can be at least partially explained by a 
series of very closely spaced frequencies around the two strongest modes. By combining 
observations taken over ^300 days we conclude that some of the closely spaced modes 
are caused by one mode whose amplitude is varying, however this frequency is still 
within '^l /iHz of an apparently stable frequency. Because of this kind of complexity 
and uncertainty we advise caution when trying to identify oscillation modes in this star. 



Key words: 

PG 1605-^072 



stars: interiors — stars: oscillations — subdwarfs — stars: individual: 



X 



1 INTRODUCTION 

The discovery of pulsations in hot subdwarf B (sdB) 
stars provides an excellent opportunity to probe sdB 
interiors using asteroseismological analysis. While it is 
generally accepted that sdBs are the field analogues of 
Extreme Horizontal Branch (EHB) stars, many ques- 
tions remain about their formation and evolution. Sub- 
dwarf B stars (like their EHB counterparts) are He- 
core burning, with 20 000 T^g <40 000K and 
5.2^1ogg^6.2 (Saffer et al. 1994). 

D'Cruz et al. (1996) showed that one possible mech- 
anism for EHB (and sdB) formation is strong mass loss 
on the Red Giant Branch. They found that the mass- 
loss efficiency required to produce EHB stars does not 
vary greatly with metallicity, which might explain the 
existence of EHB stars in both metal-rich populations 



* Based on observations made with the Danish 1.54 m telescope 
at ESO, La Silla, Chile, and on observations made with the Nordic 
Optical Telescope, operated on the island of La Palma jointly by 
Denmark, Finland, Iceland, Norway and Sweden in the Spanish 
Observatorio del Roque de Los Muchachos of the Institute de 
Astrofisica de Canarias. 



(some clusters and possibly elliptical galaxies; Dorman 
et al. 1995) and metal-poor populations (e.g., Whitney 
et al. 1994). Binary interaction is one possible mech- 
anism for such mass loss, an idea first introduced by 
Mengel et al. (1976). Indeed, at least '^60% of subdwarfs 
appear to have main sequence companions (AUard et al. 
1994; Jeffery & Pollacco 1998). 

Pulsating sdBs (also known as EC 14026 stars) 
typically have periods of 100-200 s, photometric ampli- 
tudes < 10 mmag, and oscillate in p as well as possibly 
5— modes. The pulsations are thought to be driven by 
an opacity bump due to the ionisation of heavy elements 
such as Fe at temperatures of ~ 2 x lO^K in the sdB 
envelope (Charpinet et al. 1997). 

Of all the known pulsating sdBs, PG 1605+072 has 
the most extreme properties, with the richest power 
spectrum (around 50 modes), the longest periods (up to 
~550 seconds), and the lowest surface gravity (logy ~ 
5.25) (Kilkenny et al. 1999). The low gravity implies 
that the star may have evolved off the EHB; this hy- 
pothesis is supported by modelling by Kawaler (1999). 

In Paper I (O 'Toole et al. 2000) we reported the 
detection of Doppler variations in PG 1605+072. The 
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Figure 1. Typical spectrum of PG 1605+072 obtained at the Nordic Optical Telescope. The Balmer sequence is visible from H/3 to H12, 
as well as the He i A4471 A line. The feature at ~4410 A is a bad CCD column. 



dominant frequencies detected corresponded to those 
found in photometry (Kilkenny et aL 1999), and we also 
found evidence of a wavelength dependence of the os- 
cillation amplitude of the Balmer lines. In this paper 
we present new radial velocity observations covering a 
much longer time span. Wc do not confirm the wave- 
length dependence, but wc do see evidence for very 
closely spaced frequencies, similar to the findings of 
Kilkenny et al. (1999). We also compare our results 
with those of Woolf et al. (2002), who observed this 
star with higher S/N and spectral dispersion, but over 
a much shorter time span (see Section 7). We will dis- 
cuss equivalent width measurements from our spectra 
in a future paper. 



2 OBSERVATIONS 

The bulk of the observations were made on 11 nights 
over a 16-day period in May 2000 (sec Table 1), using 
the DFOSC spectrograph on the Danish 1.54 m tele- 
scope at La Silla, Chile and the ALFOSC spectrograph 
mounted on the 2.56 m Nordic Optical Telescope on La 
Palma in the Canary Islands. To increase the frequency 
resolution of the amplitude spectrum and allow detec- 
tion of very closely spaced oscillation modes, wc also 
obtained observations at La Silla for about 1 hour per 
night in March- April 2000, before the main run (see 
Table 2). 



Table 1. Spectroscopic observations of PG1605+072. LS 
Silla; NOT = Nordic Optical Telescope. 



La 



UT-date 


Telescope 


No. of 


No. of 






hours 


spectra 


2000 May 11 


Danish 


6.87 


370 


2000 May 12 


Danish 


6.45 


300 


2000 May 16-17 


NOT 


7.31 


337 


2000 May 17-18 


NOT 


8.12 


376 


2000 May 18 


Danish 


6.75 


319 


2000 May 19 


Danish 


6.77 


375 


2000 May 19-20 


NOT 


8.40 


473 


2000 May 20 


Danish 


6.61 


370 


2000 May 21 


Danish 


2.69 


146 


2000 May 22 


Danish 


6.69 


374 


2000 May 23 


Danish 


6.67 


369 


2000 May 25 


Danish 


6.62 


346 


2000 May 26 


Danish 


1.24 


70 


Total 




81.19 


4225 



The La Silla data (March to May) consisted of 
single-order spectra projected onto a 2K LORAL CCD. 
Pixel binning (to reduce readout noise) and windowing 
(to reduce readout time) gave 66 x 500-pixel spectra 
with a total wavelength range of 3650-5000 A and a 
dispersion of 1.65 Apixcl^^. The resolution was ^6 A, 
set by a slit width of 1.5 arcsec. The exposure time was 
46 s, with a dead time of about 16 s. The average num- 
ber of photons per A in each spectrum was about 1700. 
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Table 2. Spectroscopic observations of PG1605+072 in March- 
April 2000. 



0.6 

Fractional Julian Date 



UT-date 


No. of 


No. of 




hours 


spectra 


2uUu March 18 


0.30 


12 


2uUu March zu 


1.43 


77 


2000 March 21 


0.45 


25 


2UUU March 2z 


0.45 


25 


2UUU March 26 


0.57 


32 


zUUU March 24 


0.65 


37 


2000 March 25 


0.72 


41 


20U0 March 2o 


0.93 


45 


2000 March 27 


0.93 


52 


2000 March 28 


0.83 


45 


2000 March 29 


1.17 


59 


2000 March 30 


1.05 


59 


2000 March 31 


1.07 


59 


2000 April 01 


1.15 


57 


2000 April 04 


1.03 


65 


2000 April 05 


1.42 


77 


2000 April 06 


1.47 


80 


2000 April 07 


1.55 


85 


2000 April 08 


1.65 


90 


2000 April 09 


1.57 


86 


2000 April 10 


1.70 


95 


2000 April 11 


1.90 


99 


2000 April 12 


1.23 


65 


Total 


25.22 


1367 



Figure 2. Velocity curve for H7 from May 2000. No instrumental 
variation has been corrected for. There is a 3 day gap between 
the 2nd and 3rd panels from the top. 

A similar setup was used in previous observations (see 
Paper I). 

The NOT data wore similar, with 66 x 500 binned 
pixels covering a wavelength range of 3550-5000 A. The 
dispersion was 1.65Apixel~^ and the resolution was 
'--^8 A, set by a slit width of 1.0 arcscc. The exposure 
time was 50-70 s, with a dead time of only 5 s. There 
was an average of about 8000 photons per A in each 
frame. The data from both telescopes were timcstamped 
in Modified Julian Date by the telescope computers to 
an accuracy of less than one second. 



3 REDUCTIONS 

Bias subtraction, flat fielding and background light sub- 
traction were done using IRAF, and ID spectra were ex- 
tracted using an Optimal Extraction algorithm (Horne 
1986). A cubic spline was fitted to the continuum level, 
and the spectra were normalised to a continuum value 
of unity. 

The spectrum of PG1605+072 is dominated by 
Balmer lines (see Figure 1). Also present arc some weak 
lines from He I (e.g, A4471 and A4922) and He II (e.g, 



A4686). The Balmer lines are ideal to determine Doppler 
variations. A template spectrum was created by aver- 
aging 20 high-quality spectra from each night. By using 
a different template for each night, we are effectively 
applying a high-pass filter to the data. 

To avoid pixelization effects in the cross correlation 
analysis, all spectra, including templates, were oversam- 
pled by a factor of 40 (using linear interpolation), and 
all but the Balmer line in question (one of H/?, H7, Hi5, 
He and H8) was cropped from each. 

These spectra were modified further by subtracting 
1.0 such that the continuum level was approximately 
zero. The result was passed through a half-cosine filter 
to smooth the transition to zero at the ends. 

Each resulting template was then cross-correlated 
with its respective spectra to obtain pixel displace- 
ments. The velocities were obtained by using the known 
wavelengths of the Balmer lines to convert this displace- 
ment into wavelength and then Doppler shift. See Pa- 
per 1 for further details. 

The velocity curve for H7 is shown in Figure 2. It 
has not been corrected for slow instrument drifts. The 
high quality of the NOT data (3rd, 4th and 7th panels) - 
due to its larger aperture - suggests that fainter and/or 
shorter-period pulsating sdBs could be observed with 
this telescope. 
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Figure 3. Amplitude spectra of five Balmer lines. No amplitude variation is evident in the dominant modes. 



4 TIME-SERIES ANALYSIS 

Because the quality of the observations varies through 
the data set, we performed frequency analysis using 
a weighted Fourier Transform (Kjeldsen & Frandsen 
1992), where weights were assigned to each night in 
the May 2000 data set according to the internal scat- 
ter. Weights for the March- April data were assigned to 



groups containing approximately the same number of 
observations as one night in the May data. 



4.1 Velocities of Individual Balmer lines 

The amplitude spectra of the five velocity time series 
from May 2000 are shown in Figure 3. There does not 
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appear to be any substantial variation in velocity am- 
plitude between the Balmer lines. We have analysed 
all five time scries, using the non-linear least-squares 
multi-frequency fitting software Period98 (Sperl 1998). 
The frequencies found are shown in Table 3. The same 
technique was applied to the data presented in Paper I 
(July 1999) and these frequencies are shown in Table 4. 
White noise levels in the amplitude spectra are typi- 
cally '^lOOOms"'^ for each line in the July 1999 data, 
and 600-700 ms-i in the May 2000 data. 

The major difference between the two sets of fre- 
quencies is the detection of the peaks at 1985 /iHz and 
2075 /zHz in the May 2000 data set. In Paper I it was un- 
certain whether the non-detection of the latter (found 
to be the dominant mode in photometry by Kilkenny 
et al. 1999) in any of the individual Balmer line time- 
series, was due to variable amplitude or beating between 
closely spaced modes. With a slightly longer timc-scries 
and reduced noise, this mode is clearly detected in all 
Balmer lines. Somewhat surprising is the detection of 
the mode at 1985 /iHz, as this was found in photome- 
try to have an amplitude of only 3.3 mmag, whereas the 
four other modes detected had amplitudes greater than 
13.9 mmag. This mode may have variable amplitude or 
may be a scries of unresolved modes. There is evidence 
for the latter possibility from Kilkenny et al. (1999), 
who found three modes around 1985 /iHz. As discussed 
further below, it is clear that the crowded frequency 
spectrum of PG 1605-1-072 severely complicates the in- 
terpretation of the observations. 

In Paper I we reported the possible wavelength de- 
pendence of Balmer-line velocity amplitudes (see Fig- 
ure 3 of that paper). Here, we have tried to confirm 
this result. Firstly, wc found the weighted average fre- 
quencies of each of the four modes that were detected 
in all Balmer lines. These frequencies were then simul- 
taneously fit to each time-series, giving amplitudes and 
phases for each mode. 

Figure 4 shows Balmer-line velocity amplitudes as a 
function of wavelength for these four modes. Error bars 
indicate the white noise level in the 3000-5000 /uHz re- 
gion, where there appears to be no excess power. For 
three of the modes, there is no clear evidence for any 
wavelength dependence. For the 2102 iiKz mode (tri- 
angles), all lines except H/3 appear to have the same 
amplitude, with the H/3 amplitude being slightly lower. 
We do not feel this difference is significant enough to 
claim wavelength dependence, particularly as this peak 
may in fact be a combination of several closely spaced 
modes. 



4.2 Combined Time-series 

To reduce noise and perhaps detect more frequencies, 
wc have combined the time-series for the five Balmer 
lines. Before doing this, we checked whether the oscilla- 
tion phases were the same for each mode in each time- 



4400 
Wavelength (A) 



Figure 4. Velocity amplitude as a function of wavelength. Dia- 
monds represent the mode at 2742 //Hz, triangles for 2102 //Hz, 
squares for 2075 //Hz and crosses for 1985 //Hz. Small offsets to 
some points have been added in the wavelength direction to im- 
prove clarity. There is no clear dependence of amplitude on wave- 
length. 
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Figure 5. Oscillation phase as a function of wavelength. Symbols 
are the same as in Figure 4. Small horizontal offsets have been 
added for clarity. 



series. This is simply a matter of examining the phases 
found at fixed frequency. Relative phases are plotted as 
a function of wavelength for each mode in Figure 5, with 
errors determined using simple complex arithmetic (see 
Baldry et al. 1998). There is no evidence for phase vari- 



Table 3. Frequencies (in //Hz) for five Balmer lines found from 
the May 2000 data. 



H/3 


H7 


US 


He 


H8 




1891.26 






1891.34 


1985.73 


1985.68 


1985.58 


1985.73 


1985.79 


2075.64 


2075.72 


2075.68 


2075.68 


2075.68 


2102.21 


2102.20 


2102.17 


2102.18 


2102.24 


2742.71 


2742.68 


2742.64 


2742.70 


2742.69 
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Figure 6. Amplitude spectra for July 1999 (top) and May 2000 (bottom) based on a combined Balmer line time series. Note the apparent 

change in amplitudes. 



Table 4. Frequencies (in /xHz) for five Balmer lines found from 
the July 1999 data. 

up H7 H5 He H8 

1891.02 1891.00 1890.91 
2102.14 2102.14 2102.11 2102.12 2102.12 
2742.73 2742.47 2742.69 2742.67 2742.61 



ation with wavelength, leaving us free to combine the 
time series. 

To do this, we simply averaged the velocity at each 
observation. This was possible since the weighting for 
each Balmer line was approximately the same. Weights 
were again derived from the internal scatter of each 
night. Combined time series were constructed, for both 
July 1999 and May 2000. The two resulting amplitude 
spectra are shown in Figure 6. There appears to be am- 
plitude variation over the 10 month interval between the 
observations, which we discuss in Section 5. The white 
noise level is now 555 ms^^ for the July 1999 data set 
and 425 ms~^ for the May 2000 data set, an improve- 



ment of ~45% over the best single-line amplitude spec- 
tra. 

The fitted frequencies are shown in Table 5. The 
third column contains the differences in frequency be- 
tween the two sets of data. The lengths of our two 
sets of observations are ~ 7.9 x 10^ s (July 1999) and 
~ 1.3 X 10^ s (May 2000). If the formal frequency resolu- 
tion is given by 1/T, then the corresponding resolutions 
are ~1.26/LtHz and ~0.77/LtHz, respectively. Comparing 
these values with those in the third column of Table 5, 
we see that we have been able to estimate the frequency 
to much better than the formal resolution. In the final 
column we show the frequencies found from photome- 
try, and again we see agreement that is better than the 
formal resolution. This is not surprising: at high signal- 
to-noise, the frequency of a strong peak in the power 
spectrum can be estimated to a precision many times 
better than the formal frequency resolution. 

4.3 Extended Time-series 

To improve the frequency resolution still further, we 
now add the data from March-April 2000 (Table 2) 
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Table 5. Frequencies found after combining all Balmer lines in 
the July 1999 and May 2000 data sets. Note that Av = 1^2000 — 
i^i999- The final column shows frequencies found by Kilkenny et 
al. (1999) for comparison. 



1^2000 


i^l999 


Au 


Kilkenny et al. 


{/"Hz) 


(mHz) 


(mHz) 


(1999) 


1891.19 


1890.98 


-0.21 


1891.42 


1985.75 






1985.32 


2075.68 


2075.29 


-0.39 


2075.76 


2102.21 


2102.15 


-0.06 


2101.65 


2742.71 


2742.63 


-0.08 


2742.72 




0,0 



2100 2101 2102 2103 2104 2105 

Frequency (/iHz) 

Figure 7. Spectral window (in amplitude) for 1999-2000, March- 
May 2000, March-April and May data sets in the range 2100- 
2105 fiHz. Alias peaks are seen in the top two panels with a 
~0.25/iHz spacing, and in the top panel. 

to the May 2000 series. For the extended time-series 
(T ~ 5.9 X 10^ s) the formal frequency resolution is 
~0.17/xHz. The longer time-series does add complica- 
tions however, since the number of alias peaks in the 
amplitude spectrum increases significantly. There is a 
gap of 45 days between the midpoints of the two time 
series, leading to aliasing that splits each peak into mul- 
tiplets separated by ~0.25 /xHz. The spectral window of 
each time series is shown in Figure 7. This includes the 
window function of all the 1999-2000 combined, which 
is discussed in Section 6. The extra aliasing is evident 
in the top two panels. 



The frequencies and amplitudes determined from 
the extended (March-May 2000) time-series are shown 
in Table 6. The white noise level for this time-series 
is 370 ms~^, and we use this as an error for the ve- 
locity amplitudes in column 3. The peaks at 2102 and 
2742 /xHz, both shown as single in Tables 3 and 4, are 
now resolved into three and two peaks, respectively. 
Also shown are the differences in frequency when com- 
pared with photometry. If we take into account both 
our frequency resolution and the ambiguity from alias 
peaks, the frequencies we find are in excellent agreement 
with those found in photometry (Kilkenny et al. 1999). 
In particular, the multiplicity of the peaks at 2102 and 
2742 /iHz agrees with the extended time series (T ~33 
days) used by Kilkenny et al. (1999), who found 4 fre- 
quencies around 2102 /xHz and 4 around 2742 /xHz. They 
also found 5 frequencies around 2075 /iHz, 6 around 
2270 /iHz, and 3 around 1985 /iHz. We have not de- 
tected these extra frequencies. It is not clear whether 
these multiple peaks are real or whether they caused by 
amplitude variability. We will for the moment consider 
them to be real, but will address this question in more 
detail in Section 5. 

In Table 6, we note that the peaks at 2102.48 /iHz 
and 2765.09 /iHz are matched with quite low-amplitude 
modes found in photometry (~lmmag), chosen be- 
cause they are the closest in frequency. The latter fre- 
quency can also be identified with a higher amplitude 
mode detected by Koen et al. (1998). We note also that 
2765.09 /iHz is separated from '^2742 /iHz by twice the 
one cycle per day alias frequency (i.e. 2x 11. 57 /iHz). 
The mode at 2742.47 /iHz is also identified with a low 
amplitude mode, however it was only found in the ex- 
tended time series of Kilkenny et al. (1999). Again, we 
have matched these modes purely because they are the 
closest in frequency. This appears to give further evi- 
dence for amplitude variability. 

To further quantify this variability, we have fit the 
frequencies derived from the 2000 observations to our 
1999 data. The amplitudes from the fit are shown in the 
fourth column of Table 6, next to the amplitudes from 
the 2000 observations for comparison. Errors in velocity 
are taken to be the white noise level (555 ms~^). There 
are clearly large variations in most peaks (even more so, 
considering that the weakest 1999 amplitudes should be 
considered as upper limits). 



5 SIMULATIONS 

How well resolved are the closely spaced multiplets at 
2102 /xHz and 2742 /iHz shown in Table 6? To answer 
this, we show the regions in question in Figures 8 and 
9. From the top panel of Figure 8 it is clear that the 
frequency resolution is sufficient to resolve the power at 
2102 //Hz into at least two peaks, perhaps three. The 
two frequencies found around 2742 /xHz also appear to 
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Table 6. 2000 frequencies, periods and amplitudes measured 
from the extended time-series. The fourth column shows the am- 
plitudes found using the 2000 frequency solution, the fifth column 
shows the difference between the frequencies found in this paper 

and those found in photometry (Kilkenny ot al. 1999), and the 
sixth column shows photometric amplitude rank. 



(mHz) 



P 



V2000 
(km/s) 



Viggg 
(km/s) 



(mHz) 



Rank 



1891.01 


528.82 


1.99 


3.18 


-0.41 


5 


1985.75 


503.59 


4.13 


0.63 


+0.43 


8 


2075.72 


481.76 


4.27 


2.11 


-0.04 


1 


2101.57 


475.83 


3.40 


4.97 


-0.08 


3 


2102.48 


475.63 


8.47 


11.1 


+0.04 


32 


2102.83 


475.55 


3.66 


3.05 


-0.45 


2 


2269.84 


440.56 


1.75 


0.90 


-0.27 


6 


2742.47* 


364.63 


4.45 


3.57 


-0.25 


4* 


2742.85 


364.58 


7.17 


2.87 


+0.13 


4 


2765.09 


361.65 


1.97 


0.29 


-0.20 


21 



* This mode can also be identified as one of the extra frequencies 
found in an extended time-series in Table 4 of Kilkenny et al. 
(1999). 

be resolved in Figure 9. In the bottom panel of Figure 6 
it appears that the peak around 2742 /iHz is higher than 
the one around 2102 /iHz, however Table 6 shows the 
opposite. This can be explained by examining the am- 
plitude spectra (shown in Figure 10) around 2102 /xHz 
after each frequency is removed. It is still not clear, how- 
ever, whether these frequencies are real or an artifact of 
the variation in amplitude of one frequency. 

To investigate these effects further, wc have created 
simulated time series for both our sets of observations 
(1999 and 2000). The inputs are the 10 frequencies mea- 
sured from the timc-scrics with the highest frequency 
resolution (i.e, 2000) with their corresponding ampli- 
tudes. Note that each simulation had exactly the same 
input amplitudes and frequencies, and only the phases 
were randomised. We have used two different noise lev- 
els, ^eOOms"^ and ~300ms~^, as these bracket the 
noise in our observations. The same phases were used 
at both noise levels. The basic features of the two differ- 
ent groups of spectra are the same, implying that noise 
does not play a significant role in the shape of the am- 
plitude spectrum. A selection of simulated amplitude 
spectra with noise ~600ms~^ are shown in Figure 11 
using the 2000 window function. 

The most interesting feature of the majority of 
these spectra is the variability of the amplitudes of the 
peaks at ^2102 /iHz and at ^2742 /iHz. This observed 
variability does not depend strongly on the window 
function and arises from beating between the closely 
spaced modes around these two frequencies. However 
in Table 5 there is only a small difference (<0.1/iHz) 
between frequencies measured in July 1999 and May 
2000. If there were two or more frequencies beating, wc 
would expect the peaks to move by more than this dif- 
ference, since the measured frequencies span a range of 
~1.3/iHz. 





March/April 



2102 2103 
Frequency {/^Hz) 



Figure 8. Amplitude spectra for the Combined, March/April and 
May data sets in range 2100-2105 //Hz. The dashed lines indicated 
frequencies measured from the extended time series. 



Kilkenny et al. (1999) found 2 frequencies around 
2102 /iHz in both halves of their ~15 night time se- 
ries, as well as in the full series. They also found weak 
"satellite" frequencies near almost all of the highest 
amplitude peaks. We find only one peak at 2102 /xHz 
in our May 2000 observations, however weighting the 
data may degrade the frequency resolution. The middle 
panel of Figure 8 suggests that 2 peaks may be clearly 
resolved in the March/ April 2000 observations. If wc 
assume there are 3 frequencies around 2102 /iHz, wc de- 
tect 2 frequencies in ~90% of our simulations. It seems 
that our detection of only one peak may be due simply 
to chance. Around 2742 /xHz only one frequency is al- 
ways resolved, with the 2nd frequency resolved ~50% of 
the time. This suggests a minimum of two frequencies 
around 2102 /<Hz, and possibly two around 2742 /xHz. 
Two frequencies around 2742 //Hzis supported by Reed 
(2001), who finds a similar separation (^O.SfjMz) in a 
re-analysis of the Kilkenny et al. (1999) observations. 
By combining all of our data (1999-2000), we are in a 
position to make one final test. 
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Figure 9. Same as Figure 8 expect in the range 2740-2745 /iHz. 

6 COMBINATION OF ALL 
OBSERVATIONS 

If we combine all of our observations (1999-2000), we 
find that the amplitudes of each peak is approximately 
the same (~8kms~^). This implies that the oscilla- 
tions are coherent over this one year period, and we 
can therefore perform a frequency analysis on the time 
series. Combining all of the data introduces fine struc- 
ture in the spectral window, as shown in the top panel 
of Figure 7 in the region around 2100 /xHz. The fre- 
quencies we measure from this time series are the same 
as those found in the March-May 2000 data, except 
around 2102.5 fjRz, where we now find 3 frequencies in- 
stead of two. These frequencies, while they describe the 
observations (see Figure 12), are split such that we con- 
sider them to be caused by amplitude variation. The 
two frequencies not found in 1999-2000, but found in 
the March-May 2000 data (2269.84 and 2765.09 ^Hz), 
wc consider marginal detections. We can now conclude 
that the frequencies in Table 7, derived from the March- 
May 2000 observations, are real. 

7 COMPARISON WITH WOOLF ET AL. 
(2002) 

Our observations wore by coincidence concurrent with 
those of Woolf et al. (2002), so we are in a position to 



Figure 10. Amplitude spectra of the 2100-2105 /iHz region show- 
ing the successive removal of 3 closely spaced frequencies. 



Table 7. Detected frequencies and periods in velocities of 
PG 1605+072. These are based on the March-May 2000 obser- 
vations, but also found in a time series of all observations (1999- 
2000). 



Frequency (/iHz) Period (s) 



1891.01 


528.82 


1985.75 


503.59 


2075.72 


481.76 


2101.57 


475.83 


2102.48* 


475.63 


2742.47 


364.63 


2742.85 


364.58 



*Wc include this frequency as representative of a mode with vari- 
able amplitude around 2102.5 /iHz. 

examine some of their findings. Firstly, based on data 
taken over 16.3 hours across a 32.1 hour period, they 
found an apparent power shifting between modes, which 
we also find, and attribute to both amplitude variability 
over time and beating between closely spaced frequen- 
cies. They also found evidence for rotational splitting of 
the 2742 fiRz mode, and derived a rotational period of 
12.6±2.8 hours. This is not in agreement with Heber et 
al. (1999), who derived an upper limit on the rotational 
period of 8.7 hours from spectral analysis. While the 
pulsations cause some broadening of the spectral lines. 
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Figure 12. Same for Figure 10, except with combined 1999 and 
2000 time series. 



Figure 11. Simulated amplitude spectra using the 10 frequencies 
from Table 6 as input. These spectra are based on the window 
function of all observations from 2000. It is clear that the variable 
amplitude of the highest peaks is due to beating. 



amplitude spectra. We feel timing errors may be a possi- 
ble explanation for the splitting, and lack of this central 
peak. 



the pulsation amplitudes are not high enough to shift 
Heber et al.'s upper limit into the range found by Woolf 
et al. The rotational splitting they found is ~11 /zHz, 
which is very close to one cycle per day. Because of this, 
the reality of this splitting should probably be treated 
with scepticism. 

Some of our observations overlap with Woolf et 
al.'s, and we have plotted in Figure 13 the amplitude 
spectrum of the two nights closest in time to theirs. Our 
observations for this period start at MJD 51675.1 and 
finish at MJD 51676.4, while the Woolf et al. start at 
MJD 51675.9 and finish at MJD 51677.25. The dashed 
lines in the figure indicate the frequencies found by 
Woolf et al, which appear to be one cycle per day alias 
peaks. The dotted lines indicate the frequencies found 
in this paper. For comparison wc have shown the spec- 
tral window of our observations in the top panel. The 
dashed lines correspond very well to the alias peaks 
in our spectral window, and wc conclude that this is 
most likely what Woolf et al. have identified, although 
the central peak (which we have identified with the fre- 
quency found in photometry) is not apparent in their 



8 DISCUSSION 

The model of Kawaler (1999) is the best pulsation 

model for PG 1605-1-072 produced so far, although 
it only matches 5 of the dominant frequencies from 
Kilkenny et al. (1999) (they detected up to 55 frequen- 
cies, although some may be artifacts of amplitude vari- 
ation). Kawaler 's model predicts an equatorial velocity 
of ~130kms~^ leading to a rotational period of around 
3 hours, in agreement with the upper limit found by 
Heber ct al. (1999) of P^t < 8.7h (wsim = 39kms-i). 

What arc the possible reasons for modes which 
are so closely spaced? Rotation has already been men- 
tioned, however the minimum rotational splitting for 
an / = 1 mode is ^ 16 /iHz to first order, and for an 
I = 2 mode is ^ 27 /iHz. This is based on the maximum 
rotational period of Heber ct al. (1999), and assumes 
that the star oscillates in g-modes (Kawaler 1999) in a 
similar way to a white dwarf. This is not the case, how- 
ever, since unlike white dwarfs, horizontal motion does 
not dominate vertical motion on the surface of sdBs; in 
fact, they are of the same order (Kawaler, private com- 
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Figure 13. Amplitude spectrum of niglits MJD 51675.1-51676.4 
{bottom), showing no rotational splitting of the 2742 //Hz mode. 
The frequencies found by Woolf et al. (2002) are plotted as dashed 
lines. The peaks they correspond to arc alias peaks. For compari- 
son, the frequencies found in this work are plotted as dotted lines. 
The window function for the two nights {top) is also shown. 

munication) . Despite this, the sphtting values above are 
minima, so it is clear rotation is not the only effect caus- 
ing splitting. Strong magnetic fields also cause splitting 
of oscillation modes, however little work has been done 
to measure fields of sdBs. Finally, we feel that ampli- 
tude variability over several months can explain at least 
some of the observed splittings. More detailed modelling 
of PG 1605-1-072 using the so-called second generation 
pulsating sdB models of Charpinet et al. (2001) may 
help to answer these questions. 



9 CONCLUSIONS 

Of all the pulsating sdBs currently known, 
PG 1605+072 is clearly one of the most interest- 
ing. Using our newest observations we make several 
conclusions. 

• From our main campaign we have detected 7 oscil- 
lation frequencies from velocity variations in the star. 
These frequencies are shown in Table 7. Most of these 
correspond to the strongest oscillation frequencies found 
in photometry. 

• We find the velocity amplitudes and phases in the 
individual Balmer lines to be equal and cannot con- 
firm the wavelength dependence of amplitude seen in 
Paper I. 

• Apparent amplitude variability between our two 
sets of observations can be explained by a combination 
of beating between closely spaced modes and amplitude 
variability. 

• To fully resolve and identify as many closely spaced 
frequencies as possible, it is necessary to have both good 
time coverage to reduce aliasing effects, and most im- 
portantly for PG 1605+072, to have a time series that 



is as long as possible, thus improving the frequency res- 
olution in the amplitude spectrum. 

• We can explain the possible rotational splitting of 
the 2742 /xHz peak claimed by Woolf et al. (2002) as 
an aliasing effect (possibly caused by incorrect time 
stamps), since our observations are concurrent. 

PG 1605+072 offers astcroseismologists several 
challenges. It has many oscillation frequencies, although 
some have variable amplitudes; it has longer periods 
than most other pulsating sdBs, however it is rotating 
fast enough that rotationally split modes may be un- 
equally spaced; it may be oscillating in p or g'-modes or a 
combination of both. Further theoretical work is needed 
in these areas. We feel that despite these challenges, a 
campaign of simultaneous photometry and spectroscopy 
will help to answer the questions we have raised here, 
and perhaps allow us to comprehend one of the least 
understood phases of stellar evolution. 
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